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COMPUTATIONAL METHOD 

The first-principles calculations are performed based on 
the density functional theory (DFT) [1] with generalized gra- 
dient approximation (GGA) in the form of Perdew-Burke- 
Ernzerhof (PBE) functional [2] as implemented in Vienna Ab- 
initio Simulation Package (VASP) (j. The GGA +U method 
which treats the on-site repulsion interactions of 3 J electrons 
in a mean field manner is used to evaluate the strong correla- 
tion eff'ect in TM, and a typical value ofU = 4 eV and / = 0.9 
eV are used for all TM concerned [4]. The lattice constant 
a = 3.86 A of silicene and the buckling distance d = 0.44 A 
are obtained corresponding to the global minima on the Born- 
Oppenheimer surface, which agree with existing theoretical 
data [5, 6]. As to the structure relaxation, all atoms are al- 
lowed to relax freely along all directions and all parameters 
are chosen to converge the forces to less than 0.01 eV/A. A 
vacuum space of 20 A is set to prevent the interaction between 
silicene and its periodic images along c-axis. Convergence 
tests with respect to energy cutoff' and k points sampling are 
performed to ensure numerical accuracy of total energy. We 
find an energy cutoff' of 420 eV and F centered Monkhorst- 
Pack grids of 8 x 8 x 1 for k point sampling are enough to 
converge the diff'erence in total energy to around ImeV. 



IMPACT OF STRONG CORRELATION EFFECT ON 
ADSORPTION STRUCTURE AND MAGNETIC MOMENTS 

For the sake of comparison, we define bond lengths 
{dsi-Tu) as the distance between adatom and the 3 nearest 
neighbour Si atoms to it, adsorption height ihsi-Tu) as the dis- 
tance between adatom and the lower sublattice (averaged) in 
c-axis for all 3 adsorption sites, height above Si atom (Ssi-tm) 
as the distance between adatom and the Si atom underneath 
for Ta and T^ sites (Fig. 1). 
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FIG. 1. Schematic representation of bond parameters (dsi-TM and 
^si-Tu) for 3 adsorption sites (a) H, (b) T^, and (c) T^, respectively. 



FIG. 2. The adsorption energy (A£^) of all M transition metals (TM) 
adsorbed on the three high symmetric sites of silicene monolayer, H, 
Ta and T^, in the GGA level. The inset shows the distortion energy 
{6E) of silicene when different TM are situated on the stable site H. 



In the GGA case, as can be seen from Fig. 3(a), the bond 
lengths and adsorption height generally decrease with increas- 
ing of atomic number when TM adsorbing on H site. How- 
ever, the adsorption energy doesn't follow this trend (Fig. 2), 
which has minimal value of 2.44 eV for Cr and maximal value 
of 4.75 eV for Ni. We can see from inset in Fig. 2 that the 
above trend of adsorption energy is correlated to the diff'erent 
distortion of silicene, which is energetically characterized by 
the distortion energy defined as 6E = E^is- Eg, where E^is is 
the energy of silicene after adsorption. The distortion energy 
are small for V and Cr, suggesting relatively weak interactions 
between these adatoms and silicene, therefore, the adsorption 
energy decrease from V and reaches minimal value at Cr. 

As can be seen from Fig. 5, when Sc, Ti, Cr adsorbing 
on the H site, the density of states (DOS) show peaks at the 
Fermi level, suggesting a possible Jahn-Teller distortion. In 
the case of Sc-silicene, we artificially move one of three Si 
atoms nearest to Sc to break the C3 rotational symmetry. After 
relaxation the 3 nearest Si atoms to Sc which originally copla- 
narly arranged themselves as a regular triangle (dsisc equal 
to 2.62A.) now distort to a isosceles triangle (dsiSc become 
2.63 A, 2.63 A and 3.12 A.) by pushing the moved Si atom 
down away from the upper sublattice plane by 1.24 A. The 
distorted Sc-silicene system becomes more stable than Csv 
symmetric one by lowering the total energy by O.leV. Simi- 
lar to Sc-silicene case, we could expect Jahn-teller distortion 
to further stabilise Ti-silicene and Cr-silicene systems, never- 
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FIG. 3. Structural and magnetic properties of 3 J transition metals 
adsorbed on the Hollow site of 4 x 4 silicene based on (a) GGA 
and (b) GGA+U. For Mn in GGA+U case, the bond parameters 
are corresponding to the stablest adsorption site (marked in red 
pentacle), and the magnetic moment of Hollow and sites are given 
for comparison. 



theless, the distortion for these two system turned out to be 
rather weak (The modification of position of all atoms is less 
than 0.005 A, and total energy of distorted system which no 
longer respect the Csv symmetry is lower by ~ 3meV). 

The resulting magnetic moments and magnetic instability 
aforementioned can be attributed to the interplay among crys- 
tal field splitting, the spin splitting, as well as the electron oc- 
cupation number of isolated 3d adatoms. For Sc-silicene, the 
spin splitting for Sc is relatively weak (around 0.2 eV), which 
is smaller than ligand field splitting between E2 and Ai, and 
totally there are 3 electrons occupying 3d orbitals (see Fig. 5). 
Therefore, two of these 3 electrons occupy the majority E2 or- 
bitals and the other one occupies doubly degenerate minority 
E2 orbitals, leading to I jjb magnetic moment and potential 
Jahn-Teller distortion discussed above. For Ti-silicene, ow- 
ing to the relatively large splitting of Ai (leV), the majority 
Ai orbital is occupied before the doubly degenerate minority 
E2 orbitals as indicated in inset of Ti's projected density of 
states (PDOS) from GGA (Fig. 5), resulting in peaks at Fermi 
level and magnetic moment of 2 yUg. For V, the spin splitting 
is much larger and high spin state with 5 jib moment state is 
realized, which is essential to the realization of QAHE in sil- 
icene as has been discussed before. The other cases can be 
understood in similar arguments. 

When turning on the strong correlation eff'ect, the equilib- 



rium structure of adatom- silicene systems are strongly altered 
compared with GGA case. As can be seen from Fig. 3(b), the 
dsi-TM, hsi-TM and Ssi-tm for all adsorbates (except Ni) are 
increased, especially for Ti, Cr, Mn, Fe (the bond lengths for 
these atoms increased by ~ 0.1 A while for others by ~ 0.05 
A, and the adsorption height also showed noticeable rise for 
these atoms). And the H site is still favored by most 3d TM 
(except Mn, which energetically favors T^ site). 

The change of adsorption geometry of adatom- silicene sys- 
tem can be attributed to the direct consequence of on-site 
Coulomb interactions among 3d electrons. In the case of Sc- 
silicene, the GGA + U predicts ~ 3 yUg magnetic moment 
compared with ~ 1 yUg in GGA case (Fig. 3). And the rea- 
son for this is the enhanced spin splitting, which results from 
the eff'ective U, makes spin-up A 1 and E2 states of Sc occupied 
as indicated in PDOS of Sc in Fig. 6. For Ti-silicene system, 
the net magnetic moment is ~ 3 jib rather than ~ 4 yUg as one 
may expected, indicating the 4s shell of Ti is empty and about 
one electron is transferred from Ti to silicene. This is justified 
by the lowering of Dirac cone at K and -K, though slightly 
distorted, beneath the Fermi level, and occupied PDOS of Si 
atoms closest to Ti. Similar to Ti-silicene, the Cr-silicene sys- 
tem retained ~ 5 yUg magnetic moment after transferring ^ le 
to silicene, which shifts the Dirac cone down by ~ 0.4 eV rel- 
ative to Fermi level. In the case of Mn-silicene, especially, the 
Ta site is favored if the strong correlation eff'ect is turned on, 
which is 0.07 eV (0.13 eV) lower in total energy than H (T^) 
sites. For V, Fe, Co, Ni, +U mainly enhances the spin split- 
ting while leaving the electron distribution nearly unaff'ected 
compared with GGA case, hence the same magnetic moments 
(Fig. 3). 
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FIG. 4. The evolution of band structure around valley K from the 
interplay between exchange field M and staggered potential A (in 
unit of t). The red (black) lines are for the majority (minority) spin. 
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FIG. 5. PDOS of all 3d transition metals adsorbed on the stable site (Hollow) of monolayer silicene from GGA, where positive (negative) 
values are for majority (minority) spin. The inset (b) in Sc indicates the Jahn-Teller distorted PDOS. The Fermi energy is set to eV. 
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FIG. 6. PDOS of all M transition metals adsorbed on the Hollow (H) site of monolayer silicene from GGK+U (with an effective U being 
S.l^V), where positive (negative) values are for majority (minority) spin. For Mn, the PDOS including 4^ orbitals are given in inset when 
absorbing on the stable site T^. The Fermi energy is set to eV. 



BAND STRUCTURE EVOLUTION FROM TIGHT BINDING 
MODEL 



Fig. 4 illustrates the evolution of the band structure around 
valley K with different combination of exchange field M and 
staggered potential A. Panel(a) shows the band structure of 
pristine silicene with perfect Dirac-like energy dispersion. In 
panel(b), the spin degeneracy is lifted when only exchange 
field M is turned on, rendering crossing of bands with op- 
posite spin around K point in the reciprocal space similar to 
graphene case |7]. While if only staggered potential A is ap- 



plied (panel(c)), the system becomes insulating with the va- 
lence and conduction bands twofold degenerated. 

As mentioned before, in ?>d TM doped silicene, the ex- 
change field and staggered potential should appear simultane- 
ously. Panels(d)-(e) indicate the interplay between exchange 
field (M) and staggered potential (A). In the case of panel(d) 
where M = A, there always exists a degenerate point right 
at the Fermi level. When increasing M and making M > A 
(panel(e)), the two spin subbands near Fermi level cross, re- 
sulting a circular Fermi surface. While if keeps increase ex- 
change field M (panel(f)), the system enters insulating state. 
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The M > A case is of particular interest in our investigation, 
since upon turning on the spin orbit coupUng effect, the case 
(e) will give rise to the QAHE and electrically tunable topo- 
logical states from QAHE to QVHE as discussed before. 
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Introduction -The recently discovered topological insula- 
tors (TIs) m-El] have aroused great interest in the fields of 
condensed matter physics and materials science due to the 
metallic boundary states protected by time-reversal symme- 
try (TRS). TIs have also become perfect breeding ground for 
a variety of exotic quantum phenomena In particu- 

lar, breaking the TRS respected by TIs via magnetic doping 
is predicted to give birth to the Majorana fermion [6J, topo- 
logical magnetoelectric eflTect JtI], and the so-called quantum 
anomalous Hall effect (QAHE)r8'-n^, which has quantised 
Hall conductance in the absence of external field [13] and can 
be intuitively thought as half of TIs. Since TIs have been fabri- 
cated in materials ranging from 2D [14] to 3D [l?], engineer- 
ing these novel phenomena in real materials represents one of 
the most fascinating areas in this field. 

As the 2D TI (also known as quantum spin Hall (QSH) 
insulator) graphene Silk] has been altering the research 
direction of nanoelectronics from silicon-based materials to 
carbon-based ones, however, the advent of silicene [UtI - EiI] . 
which is the silicon equivalent for graphene, seems to turn the 
tide. Silicene is closely analogous to graphene in the sense 
that it consists of a single layer of Si atoms arranged to a 
low buckled honeycomb lattice, and its low energy physics 
can be described by Dirac-type energy-momentum dispersion 
akin to that in graphene [17], hence the inherited many in- 
triguing properties, including the expected Dirac fermion s and 
QSH effect 11221] . Yet a striking difference between silicene and 
graphene is that the stable silicene monolayer has additional 
buckling degree [17], which accounts for the relatively large 
(1.55meV) spin orbit coupling (SOC) induced gap 1E2I] in sil- 
icene and a couple of unusual quantum phenomena recently 
reported [23-25]. Indeed, these features together with the nat- 
ural compatibility with current silicon-based microelectronics 
industry make silicene a promising candidate for future nano- 
electronics application. Moreover, from the view of practical 
applications, it is highly appreciated if magnetism or sizable 
band gap or both, like in QAHE with additional edge states 
protected by topology, can be established in the nonmagnetic 
silicene, especially in the presence of the bulking degree. 

In this Letter, we present a systematic investigation of 
the adsorption properties and magnetism of 3d transition 



metals (TM) doped silicene. We demonstrate that 3d TM 
strongly bonding with silicene and the TM- silicene systems 
are strongly magnetic. From combined tight-binging (TB) 
model analysis and first-principles Wannier interpolation, we 
investigate the topological properties of the resulting TM- 
silicene systems. Our results suggest the V doped silicene 
hosts a stable QAHE, and this system can also be half-metallic 
[26] if the Fermi level is properly tuned. Further, a close study 
of the TB model in the band inverted regime gives rise to 
another topologically nontrivial state, which supports quan- 
tum valley Hall effect (QVHE) We predict the result- 
ing QAHE and QVHE can be tuned directly using an external 
electrical field. 

Adsorption and magnetism analysis .- We use 4x4 super- 
cell of silicene to model the interaction between 3d TM (Sc, 
Ti, V, Cr, Mn, Fe, Co, Ni) and silicene. As silicene has buck- 
led geometry, we consider three high symmetry adsorption 
sites, namely the hollow (H) site at the center of a hexagon, 
two top sites denoted as T^ and Tg corresponding to the top 
of Si atoms belonging to A and B sublattice, respectively (see 
Fig. 1). To evaluate the eflTect of on-site Coulomb interactions 
among 3d electrons of adatoms on the equilibrium structure 
and magnetic properties of the TM-silicene system, the simu- 
lations have been carried out within generalized gradient ap- 
proximation (GGA) [31] and GGA-\-U [32] framework sepa- 
rately. 

Let us first focus on the GGA case. From our first- 
principles results, of the three adsorption sites concerned, all 
3d TM energetically favor H site, which is 0.02 eV ~ 0.60 eV 
and 0.20 eV~ 0.80 eV higher in adsorption energy (A^) iH 
than Ta and Tb site, respectively (see Fig. 2 in Ref. [29]). The 
bondings between 3d TM and silicene are strongly covalent as 
manifested by a much larger AE ranging from 2.44 eV to 4.75 
eV. The unusual large AE compared with that in graphene case 
llUllll] could be related to the covalently more active sp^ -like 
orbitals of silicene, which result from the unique buckled ge- 
ometry . 

Much like the graphene case fll','33], most of TM (except 
Ni) doped silicene exhibit magnetism with sizable magnetic 
moments ranging from ~lyUfito~5yUfi. A relatively large 
magnetic moment is key to the realization of QAHE in silicene 
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FIG. 1. The lattice geometry of 3d transition metals doped silicene 
with lattice constant \a^\ - |^4| = 4a, where a - 3.86 A is the lattice 
constant of silicene. (a) Top view of 4 X 4 silicene monolayer where 
the 3 adsorption sites (Hollow (H), top of A sublattice (T^) and B 
sublattice (T^)) are marked out with black letters, (b) Side view of 
silicene, the two equivalent Si sublattices are labeled as A and B, 
respectively, with a buckled distance 6 - 0.44 A. 



which we will discuss later. We also note that when some TM 
(Sc, Ti, Cr) adsorbing on H site, the density of states (DOS) 
show peaks at the Fermi level, indicating that these systems 
could be magnetic instable and may undergo Jahn-Teller dis- 
tortion to lower total energy (especially for Sc, see Ref. [|29ll 
for detailed discussion). 

The resulting magnetic moments and possible Jahn-Teller 
distortion aforementioned can be understood in the light of 
symmetry considerations. When TM are deposited on high 
symmetry sites of silicene (H, T^, Tg), the 3 J subshell of 
adatom split into three groups under the symmetric crystal 
field of system: the 3 ^3^2.^2 state corresponding to Ai symme- 
try group, the twofold degenerate Ei group consisting of the 
?)dxz and ?>dyz states and the E2 group consisting of the ?>dxy 
and ?>dx2-y2 states. Therefore, the three groups of 3 J states 
hybridise with tt orbitals of silicene weakly or strongly ac- 
cording to the diff'erent symmetrical properties in similar way 
as in Benzene |[34ll and graphene cases |35]. Since 3d or- 
bitals are anti-bonding states, the energy order of them are 
usually 6(^2) < ^(^1) < ^{Ei) and in general 6(^2) and e{Ai) 
are close to each other due to similar hybridization strength 
with 71 orbitals. After incorporating spin polarization, the 3 
groups of states split, according to diff'erent splitting energy, 
into 10 spin polarized orbitals. Meanwhile, the outer 4^ elec- 
trons of adatoms experience relatively large electrostatic inter- 
action from the tt manifold of silicene than 3 J shells due to its 
spherical symmetry and delocalized nature, making possible 
charge transfer from the Asio?>d shells. Thus, the crystal field 
splitting, spin splitting, together with the occupation number, 
mostly dominate the electronic structure of adsorbed TM ions 
(see Ref. fl^ for detailed discussion). 

Since the strong correlation eff'ect of 3J electrons is not 
negligible for a practical description of adsorption, we next 
consider GGA+U case. After full relaxation, the adsorption 
geometry of adatom- silicene system is strongly altered com- 



pared with GGA case (see Fig. 3(b) in Ref. Q). And 
most 3 J TM still favor H site (except Mn, which now en- 
ergetically favors Ta site). Clearly, the geometry change of 
adatom- silicene systems are the direct consequence of on-site 
Coulomb interactions among ?>d electrons, which modify the 
electron distribution in 3 J and 4^ shells of adatoms and n or- 
bitals of silicene as can be seen from the changes of magnetic 
moments (see Ref. [291 for detailed discussion). 

Effective model and Chern number analysis -In this part, 
we turn to the main finding of this work, namely the QAHE 
in the absence of external field via doping ?>d TM and the pre- 
diction of electrically tunable topological states. As has been 
shown in Refs. 1 1 1, 12], the QAHE could be realized via dop- 
ing certain ?>d or 5d adatoms on the hollow site of graphene. In 
Fe doped graphene case 1 11], the QAHE gaps occurs around 
the Dirac K points of the Brillouin zone, and the low energy 
physics can be described by a Hamiltonian for graphene in the 
presence of extrinsic Rashba SOC (/l^^O and exchange field 
(M) [36] introduced solely by the adatom. Here in silicene, 
however, owing to the low buckled structure, there exists in- 
trinsic Rashba SOC (4«0 iS. Moreover, when depositing 
3 J TM on the stable adsorption site, the induced inequality 
of AB sublattice potential (A) necessarily arise and compete 
with magnetization. Below we identify conditions for the real- 
ization of QAHE in silicene based on a eff'ective Hamiltonian 
|[3Q] by introducing a staggered AB sublattice potential be- 
sides SOC {A^^^ and A so) and exchange field (M). In the basis 
of {A, B} <S) IT, vl}, the Hamiltonian reads: 

with = SeffTo (To ± T3 /Zii flVrikxTi =F kyT2) CTq, 

= A^^\±Ti (T2 - T2 (Ti) + At3 (Tq + Mtq (T3, 

hn = -AsoO's - aA!^\kyCri - kxcr2), where H^^^ are the total 
Hamiltonian for the two inequivalent Dirac points K(-\-) and 
-K(-), Hf are the low energy eff'ective Hamiltonian for the 
QSH insulator silicene, include all the eff'ects introduced 
by the 3d dopants, including eff'ective spin-dependent mag- 
netic field M, site-dependent staggered potential A and the re- 
sulting extrinsic Rashba SOC A^j^K The r and cr are the Pauli 
matrices acting separately on pseudospin (sublattice) and spin 
space, Seff stands for the si - A2nd term in Ref. ['22'], vf and 
a are the Fermi velocity and the lattice constant, respectively, 
and Aso is the eff'ective SOC. 

By diagonalizing the above Hamiltonian, the electronic 
structure around each valley (K and -K) in Brillouin zone can 
be obtained (the Se/f term can be safely ignored). From the 
interplay between exchange field (M) and staggered potential 
(A) (see Fig. 4 in Ref. 11291] ). we conclude that energy bands 
with opposite spin intersect, which is essential to the QAHE 
in Fe doped graphene, only when M/A > 1. This is diff'erent 
from that in Ref. jj^. where A = and there always exists 2 
degenerate points around Dirac point as long as M 9^ 0. 

We found that either the extrinsic or the intrinsic Rashba 
SOC would lead to insulating state when starting from the 
case M/A > 1. To identify the topological properties of the 
resulting insulating state, we resort to the Chern number (C) 
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(a) Total Chern number (b) Chern number at K (c) Chern number at -K 
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FIG. 2. The transition of Chern number by tuning A^j^^ and A^^^ (in unit 
of 0- (a) Three topological nontrivial states, QAHE(2), QVHE(O) 
and QAHE(-2) with Chern number +2, 0, and -2, can be obtained 
from different combination of A^^^ and /l^^ (b) and (c) represent the 
the variation of Ck and C-k- (d), (e) and (f) depict the band structure 
along ky = line in BZ for the three topological states (QHE(2), 
QVHE(O) and QHE(-2)) in (a). The yellow integers (+1) repre- 
sent Ck and C-k, corresponding to the sum of topological charge 
of each valence bands (the white +0.5 and +1.5). (g) and (h) show 
the gap closing around K and -K. They are the transition states from 
QAHE(+2) to QVHE(O) and from QVHE(O) to QAHE(-2), respec- 
tively. 



analysis 1371] . The C can be obtained by the integral over the 
first Brillouin zone (BZ): C = ^ ^(k)£k. The Q(k) is the 
usual Berry curvature of all occupied states [|38|: 

<<A«k|Vxl<Amk)<<Amk|V3;|<A«k) 



(2) 



where is the Fermi-Dirac distribution function for band n, 
i//nk is the Bloch function of eigenenergy 6„k, Vx, Vy are the ve- 
locity operators. And anomalous Hall conductivity is readily 
given by cr^y = (e^/h)C. Interestingly, the extrinsic Rashba 
SOC (/l^^O gives insulating state with C = +2 while the in- 
trinsic one (/l^O leads to that with C = -2. One may expect 
that diff'erent Chern number state can be realized if tuning the 
two types of Rashba SOC properly in experiments. Fig. 2(a) 
indicates that this is indeed the case, where C can take integer 
value of +2, 0, -2 with diff'erent combination of A^^^ and A^j^K 
Careful study shows that the tunable C originates from diff'er- 
ent response of bulk gap to the 2 types of Rashba SOC (A^^^ 
and /l^^O around K and -K. When increasing A^j^^ while keep- 
ing /l^^ fixed, for example, at 0.03^ (t = 1.6 eV is the nearest 
neighbor hopping parameter [30J) as shown in Fig. 2(a), we 
can see clearly the transition of Chern number of each valley 



(Ck and C-k) from to -1 but with diff'erent rate, i.e., Ck 
experiences a topological transition earlier than C-k (see Figs. 
2(b)-2(c)). The step change of Ck and C-k is justified by the 
observation of bands touching and gap reopening around each 
valley (Figs. 2(d)-2(f)). Notice that the rotational symmetry 
of the eff'ective Hamiltonians along z direction in any angle is 
broken after bringing in the Rashba SOC terms, and hence the 
band touching happens only on = line in BZ for valley K 
(Fig. 2(g)) while kx = line for valley -K (Fig. 2(h)). Conse- 
quently, the system can be in QAHE phase (with C being -h2 or 
-2) or QVHE (with C being and Ck = -C-k = -1) depend- 
ing on diff'erent value of A^j^\ which is controllable through an 
external gate voltage. 

The eff'ective SOC (Ago), however, further breaks the 
particle-hole symmetry of the above Hamiltonian, making en- 
ergy bands shift up (at valley K) or down (at valley -K) rel- 
ative to the Fermi level while leaving the topological charge 
of each valley unchanged. Hence, as long as the shifting is 
small, the system is still insulating and the above discussion 
of topological transition remain valid. 

Some of the topologically non-trivial phases can be realized 
in 3d TM doped silicene as predicted by our first-principles 
calculations. For example, we notice that the opposite spin 
subbands cross around the Fermi level in V, Cr, Mn doped 
silicene due to relatively large magnetization. Meanwhile, for 
these systems the spin-up and spin-down subbands near Fermi 
level are also gapped (Fig. 3(al) and Fig. 3(a3)), which makes 
these systems candidates for half-metallic materials if tuning 
Fermi level properly. We take V- silicene as a prototype and 
discuss in detail, as this system is insulating with an energy 
gap around 6meV (Fig. 3(a2)) when only the SOC is turned 
on. We find the SOC induced band gap in V- silicene is stable 
even strong correlation eff'ect of V is considered by including 
eff'ective U value ranging from to 6 eV (Fig. 3(a4)). 

To obtain an accurate Berry curvature distribution and 
hence Chern number in first-principles level in V doped sil- 
icene, the Wannier interpolation technique have been used 
. The band structure from first principle is well repro- 
duced by interpolation as indicated in Fig. 3(a4). In Fig. 3(b), 
the Berry curvature (^z(k)) in /:-space is explicitly shown. As 
we may observe, the most nonzero values of Berry curvature 
distribute around the Dirac K points by forming small circles, 
where exactly the avoided crossing happens (Fig. 3(a4)). By 
the integral over BZ, we indeed find the Chern number of all 
occupied bands equals to an integer value of +2, which signals 
the V doped silicene is in QAHE phase. 

For the proposed QVHE state in former model analysis, we 
note that the A"^^ is rather small (being about 4.4 x lO'^^O 
Therefore, the QVHE region in Fig. 2(a) would be quite small 
and it may not be easy to directly observe the QVHE from 
experiment. To tide it over, putting V-silicene system or its 
counterpart of TM doped Germanium on substrate may be a 
feasible solution. 

Conclusions -To summarize, we have demonstrated that 
the 3d TM doped silicene can be intriguing materials as man- 
ifested by the induced strong magnetic moments, potential 
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FIG. 3. (a) The band structures of V doped silicene from GGA ((al)-(a2)) and GGA+U ((a3)-(a4)), respectively. The red (black) color in (al) 
and (a3) correspond to majority spin (minority spin) subbands. After including SOC effect, a gap is opened at the Fermi level ((a2) and (a4)). 
In (a4), the band structure from Wannier interpolation is also shown in pink dashed lines, (b) The distribution of Berry curvature (^^(Q) in V 
doped silicene from GGA+[/+SOC. The BZ is marked out with black hexagon. 



half-metallic property as well as sizable topologically non- 
trivial gaps. These features have also been confirmed in the 
presence of strong correlation eff'ect of 3d TM. Moreover, we 
predicted the emerging of electrically controllable topological 
states (QAHE and QVHE phases as characterized by diff'er- 
ent Chem number) in certain TM-silicene systems where the 
energy bands are being inverted. Our work may provide new 
candidate for the long- sought QAHE and platforms to manip- 
ulate topological phase transition electrically. 
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